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ABSTRACT 

While most of the microwave background fluctuations on angular scales 
greater than a few arcminutes were generated at z > 800, the low redshift 
universe does distort the microwave background. Since the Sloan Digital 
Sky Survey (SDSS) traces the structures in the low-redshift universe, we 
can gain additional insights into the physics of the low-redshift universe by 
cross-correlating microwave background maps with template maps produced 
from the SDSS. We present a formalism for cross-correlating data from the 
Microwave Anisotropy Probe (MAP) with the Sloan Survey for the thermal 
Sunyaev-Zel'dovich (SZ) effect, the Integrated Sachs- Wolfe (ISW) effect, and 
weak lensing. This formalism is used to compute the signal-to-noise for 
cross-correlating these effects with various classes of tracer objects from the 
SDSS. It is found that the anticipated samples of SDSS quasars and galaxies with 
photometrically-determined redshifts would be good tracers for cross-correlating 
with the CMB. We find that the SZ-galaxy cross-correlation would give good 
constraints on pressure fluctuations in supercluster-scale gas. Cross-correlating 
weakly-lensed quasars with maps of the convergence of the CMB is found to 
give strong constraints on Qq as well as the equation of state w. We find that 
the ISW cross-correlation gives poor signal-to-noise using these techniques. 
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1. INTRODUCTION 



Microwave background observations and measurements of large scale structure are 
complementary. While the Sloan Digital Sky Survey (SDSS; punn 1995[ ) measures the 
distribution of galaxies in the nearby universe, NASA's Microwave Anisotropy Probe 
(MAP; [Bennett et al. 1997|) primarily probes the distribution of matter and radiation at 
z=1300. Since the nearby universe is thought to have the same statistical properties as 
the distant universe, we can gain insights into low-redshift physics as well as constrain 
cosmological parameters by cross-correlating the final conditions measured by SDSS with 
the initial conditions probed by MAP ( [Boughn, Crittenden, fc Turok 1998| ; |Eisenstein, Hu 
fc Tegmark i99l ; |Wang, Spergel, fc Strauss 1999Q . 

To maximize the signal-to-noise from the cross-correlation within the limitations on the 
CMB data set by the specifications of the MAP Project, a good tracer of the low-redshift 
density distribution must fulfill two criteria: 



1. The tracer objects must be numerous enough that the Poisson error term in the 
cross-correlation is small. 

2. The tracer must probe as large a volume of redshift-space as possible so that the 
cross-correlation signal is maximized. 



In practice, one finds that object types with high number densities (e.g. galaxies) tend 
to have redshift distributions that only probe a small part of the volume and object types 
that probe a larger fraction of the volume (e.g. quasars) are less numerous. Therefore 
it is necessary to strike a balance between these two criteria in order to optimize the 
signal-to-noise. 

In this paper, we investigate several classes of objects that will be available in the 
SDSS survey in order to find tracers that optimize the cross-correlation with the CMB. 
We find that SDSS quasars with magnitude i' < 21 and galaxies with magnitude r' < 21, 
which will have photometrically-determined redshifts, are good tracers. We then use the 
cross-correlation formalism to show that, 



1. cross-correlating the SDSS galaxy sample with r' < 21 in various redshift "slices" 
with maps of the Sunyaev-Zel'dovich effect from MAP gives good constraints on the 
pressure fluctuations in supercluster-scale gas, and 

2. cross correlating the SDSS quasar sample with i' < 21 with convergence maps from 
MAP gives strong constraints on the density parameter and the equation of state. 
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The second-order effects that alter the microwave background and the low-redshift 
tracers used in the cross-correlation are introduced in § |2|. In § ^ we outline the cross- 
correlation formalism, derive an optimal filter for maximizing the signal-to-noise from it 
and discuss briefly the computational issues involved in the cross-correlation. The results of 
the calculations are presented in § |j. We explore the possibility of constraining the equation 
of state with these methods in § [5[ Finally, our conclusions are summarized in § || 

2. THE LOW-REDSHIFT UNIVERSE AND THE MICROWAVE 

BACKGROUND 

There are three primary low redshift effects that alter the cosmic microwave background 
(CMB). 

2.1. Sunyaev-Zel'dovich (SZ) Effect 

Hot gas in clusters and filaments produce new distortions in the microwave background. 
The SZ effect QSunyaev fc Zel'dovich 1972, 1980|) is the change in energy experienced by 
CMB photons when they scatter from the hot gas. It takes two forms. The dominant 
contribution is the thermal SZ effect, the gain in energy acquired from the thermal motion 
of the gas. The kinetic SZ effect arises from the Doppler shift due to the bulk motion of the 
gas. 

2.2. Integrated Sachs- Wolfe (ISW) Effect 

In low-fio universes, the fact that potentials decay (in other words, matter fluctuations 
stop growing) at some epoch makes another contribution to the large-scale CMB anisotropy. 
In addition to the redshift experienced while climbing out of potential wells on the last 
scattering surface, photons experience a cumulative change in energy due to the decaying 
potentials as they travel to the observer. The blueshift of a photon falling into a decaying 
potential well is not entirely cancelled by the redshift as it climbs out. This leads to a net 
energy change, which accumulates along the photon path. This Integrated Sachs- Wolfe 
effect is distinct from the more commonly considered redshifting at the last scattering 
surface which has become known as the Sachs- Wolfe effect. Both effects are considered in 
the original paper by |Sachs fc Wolfe (1967)| . 

The ISW effect operates most strongly on scales where the change of potential is large 
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over a wavelength. For ACDM models, the effect is confined to the largest angular scales 
( [Kofman fc Starobinsky 1985| ), i.e. I < 10's. 



2.3. Weak Lensing 

Light propagating through an inhomogeneous universe is distorted through lensing by 
foreground matter ( |Gunn 1967| ). Dark matter distributed along the line of sight between 



the surface of last scattering and the present epoch deflects CMB photons and imprints a 
characteristic pattern onto the CMB anisotropies created by acoustic oscillations. On large 
and intermediate scales lensing smoothes the acoustic oscillations while on very small scales 
it creates additional power ( ]5eljak 1996|) . 

Zaldarriaga & Seljak (1999)' have developed a useful formalism for computing the 
convergence of the microwave sky. The convergence is a quantity which induces random 
deflections in the direction of the CMB photons. It is the integrated mass density along the 
line of sight divided by the critical density. 



2.4. Low-redshift Tracers 

The SDSS, which will cover about a quarter of the sky, will contain many classes of 
objects which trace large scale structure. These can be cross-correlated with microwave 
background maps produced by MAP. Some of these are: 

1. Photometric-z galaxies (r' < 21) 

2. All galaxies (r' < 23) 

3. Spectrometric-z quasars {%' < 19) 

4. Photometric-z quasars (i' < 21) 

5. FIRST objects 

Here, the magnitude cuts are made in the SDSS r' and i' photometric bands. The 
first and fourth categories will have photometrically determined redshifts. The magnitude 
cuts for these categories reflect conservative estimates of the limiting magnitudes to 
which photo-z techniques can be applied (D. Eisenstein 2000, private communication, 
G. Richards 2000, private communication). Upper estimates for the expected errors in 
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photometrically-determined redshifts for r' < 21 galaxies are about 0.04-0.08. Typically, at 
present, photometric methods for quasars give 50% of redshifts correct to within ±0.1 and 
70-75% correct to within ±0.2. This is fairly independent of magnitude down to i' = 20. 
After that, photo-z cannot be done reliably on z < 2.2 QSOs. However, for z > 2.2 QSOs, 
the redshift errors should be comparable. These errors should improve in the near future as 
details of the photo-z techniques are refined. 

The second category contains galaxies brighter than the 5<r detection limit of the 
SDSS r' filter. Since we will not have redshifts for all of these galaxies, the results of 
cross-correlating this sample with low redshift effects are of academic interest only. The 
third category contains quasars for which redshifts would be obtained by the SDSS 
spectroscopic survey. For the SDSS spectroscopic quasar sample, the redshift errors are 
negligible (< 0.1%). The last category contains objects from the FIRST survey QWhite et 
aL 1997|) which would be optically matched by the SDSS. 

Approximating the radial distribution function of the tracer sample by a 5-function at 
the mean redshift of the sample changes cross-correlation results by 5-10%. We estimate 
the signal-to-noise to change by roughly the same factor due to photometric redshift errors. 



3. CROSS-CORRELATION FORMALISM 

In this section we introduce the notation and develop the formalism for cross-correlating 
each of the tracer object types with the secondary CMB anisotropies described above. 

We begin with the comoving Friedman equation for a spatially flat universe with a 
cosmological constant A: 



8ttG a 2 



Po , 
-3 + Pa 



a 2 H 2 



-4 + a 



n ) 



where a is the dimensionless scale factor, prime denotes derivatives with respect to 
conformal time r, Hq is the Hubble constant and G is the gravitational constant. We take 
Qo + Q\ = 1, where f2 = 8nGpo/ (3Hq) and = 8tiGpa/ (3Hq), where po and pa are the 
mean densities in dark matter and vacuum energy today, respectively. We normalize the 
scale factor so that a(r ) = 1, where r is the present conformal time. The linear growth 
factor ( [Peacock 1999, p468|) 
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is normalized so that D(t ) = 1. We now introduce the conformal lookback time 
i] = (tq — t), since it is more natural to frame the integrals that follow in terms 77. We have 
used the conformal transformation dr = (c/a)dt, where t is the standard time coordinate, 
with the result that the conformal lookback "time" represents the proper motion distance. 

The fluctuations in projected surface density of tracer objects i in a redshift slice in 
direction n ( [Moessner, Jain, fc Villumsen 1998]) are given by, 

£<(n)=/ drj[b i 6(x)W(rj)+2(2.5s i -l)K] (3) 

In this expression, the first term in the integrand gives the fluctuations on the sky resulting 
from intrinsic clustering of the tracer objects. The second term gives the fluctuations arising 
due to magnification bias in the weak lensing limit, k < 1: lensing increases the area of a 
given patch on the sky, thus diluting the number density. On the other hand, tracer objects 
too faint to be included in a sample of a given limiting magnitude are brightened as a result 
of lensing and may therefore be included in the sample. The net magnification bias can 
lead to either enhancement or suppression of the observed number counts, depending on 
the slope of the number-magnitude relation. 

Here, we have assumed that the tracer density perturbations are biased linearly with 
respect to dark matter density perturbations, i.e. <5j(x) = 6j5(x), where for simplicity, the 
bias factor bi is assumed constant and independent of scale and redshift. In reality, the 
bias parameter is likely to be time- and scale-dependent, and could be estimated for real 
galaxy/quasar samples using their power spectra and factored into the cross-correlation. 

W(rj) denotes the radial distribution of the tracer objects, is the logarithmic slope 
of the number counts of tracers N (m) in a sample with limiting magnitude m, 

_ d\ogN (m) 
dm 

k is the convergence, a weighted projection of the density field along the line of sight: 

k = -Q [ — / dr]- -<5(x) (5) 

2 \ c / Jo a[To — rj) 

where the radial window over the dark matter fluctuations 5 is g/a. In a flat Q — 1 universe, 
the radial weight function is given by, 

g( V ) = r] r di^^Wii) (6) 

Jr, T}' 

where W(r]') denotes the normalized radial distribution of the lensing objects, i.e. the 
probability of finding an object at a given distance. 
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Expanding the comoving matter density fluctuations as, 



d>(x) = S(hrj) 

d 3 k 



(2 f )3/2 1) ( r 0-#( k ) e 



(7) 



Eq. |3| becomes 
Si(n) = 



TO 



drj 



/V 2 



Mr(//) + :«>„(-) (2.5 Si -l) / (7?) 

c / a(r - 77) 



rf 3 k 

(2tt) 3 / 2 



L>(t - V )6{k)e 



We will cross-correlate projected surface density fluctuations of the quasar and galaxy 
samples from the SDSS, described in § 0, with the CMB. 

In the case of the optically matched FIRST sample, we will cross-correlate instead 
the ellipticity e = 2(1 — <j 2 )k where the RMS ellipticity error is taken to be cr e = 0.4 
( JRcfregier et al. 1998|) . The weak lensing shear 7, is related to the source-averaged 
ellipticity by (e*) ~ — gji, where g = 2(1 — a 2 ) is the shear-ellipticity conversion factor, 
and o\ = (e 2 ) = (e 2 ) is the variance of the intrinsic source ellipticities. Since we have 
71 = Kcosa and 72 = kcos/3, (•y 2 ) = (7^) + (7I) = (k 2 ) and the above expression follows. 



Thus, 



e(n) = 1.68 x ^Q Q 



T0 , g(v) f d 3 k 

(, V a(r -r l ) J (2vr) 3 / 2 



D(r - V )5(k)e^ 



(9) 



Now we perform the cross-correlation with the CMB maps for each of the three 
low-redshift effects described in 8 0. 



3.1. SZ Effect 



The thermal SZ effect (§ |2.1|) in a given direction is computed as a line-integral. In the 
Rayleigh- Jeans portion of the spectrum (i.e. the long-wavelength limit) this is given by, 



AT 



T 



m 



sz 



-2a T k B b gas r drfne{rf)Urf)m 
m e c A Jo 



-2a T k B n°T% 



TO 

dig a(ro — rj 



r\-P 



d 3 k' 

(2vr) 3 / 2 



D{r - V ')5(k')e lk ' ^' (10) 



where n° e and T e ° denote the electron number density and temperature at the present epoch, 
respectively, and (3 is a parameter describing the redshift evolution of n e and T e . For the 
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calculations in this paper, we assume that (3 = 0, i.e. that the pressure is independent of 
redshift. Note that a(r —7)') = 1/(1 + z). b gas is the bias factor of the gas relative to 
the dark matter density perturbations <5(x'). Again, we have assumed a linear bias model 
with a constant bias factor for simplicity. In order to produce a conservative estimate of 
signal-to-noise for the SZ correlation, we have used b gas = 4.0, which is somewhat lower 



than the average bias predicted by numerical simulations for the low-fc regime QRefregier et 
al. 1999|) . 



n° e is calculated assuming VL^h 2 = 0.015 where h is the Hubble constant in units of 100 
kms" 1 Mpc" 1 . T e ° is taken to be 0.5 keV. a T = 6.65 x 10~ 25 cm 2 is the Thompson scattering 
cross-section, c is the speed of light, and m e is the electron rest mass. 



Using Eqs. ||] and |n], the CMB-number count cross-correlation is 
2a T k B n Xb ga , 



AT 

— (m)Ei(n 
J- / sz 



Here, we have used 



m P c z 



drj 



iV 2 



<>M'iij) ■ 3<> n ( — ] (2.5s, -1) ( 9{r}) , 
c / a(r - 7]) 



x n^-^r^E^l^k)!^^^ (11) 

•'0 a(r - 7]') p J (27r J J 



(12) 



(5(k)8(k>)) = 5 Dirac (k + k r )P(k) 
where P(k) is the mass power spectrum. Expanding out the exponentials yields, 

J dQk^e-^' = 4tt ]T(2£ + l)j e (krj)j e (krj')P e (ih ■ n) (13) 



Thus, we have 



where 



AT 

— 1 



E 



(21 + 1) ^Sj-sz p / * ~x 
— — — i^(m ■ nj 



with 



u> £ *(/c) = / dr/ 

JO 



C, 



TO 



Si _ 5Z 2 (2a T k B n°Xb : 



71 



H^ 2 



""- 1 ' k 2 dkP{k)wf{k)w^{k) 



(14) 
(15) 



and 



W(^)+3fio — (2.5s, - 1, 

c / a(r - ?/) 



'o a(r - rf) 
Note also that the number count autocorrelation is 

Cf ~ s < = - / k 2 dkP(k) \wf{k)] 2 

71 J L J 



D(T -v)je(kr)) (16) 



;i7) 



(18) 
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The signal-to-noise per £-mode for a cross-correlation, assuming the correlation is 
weak, is (S/N)j = (ab) 2 / ({a 2 ) (b 2 )) . Applying this formula to the SZ-number count 
cross-correlation yields, 

I) 2 =*>V(2t+l) (C *~">' (19) 

where Cj T is the total CMB signal (which is dominated by contributions from the surface 
of last scatter), AQ = n is the SDSS sky coverage, and the shot noise contribution per 
mode is, 

c POISSON = 1 (2Q) 



where iVj is the number of tracer objects i per steradian (Peebles 1980, Eq. 46 ■ 1 6 | ) . The 
detector contribution per ^-mode is, 

C? ET = nKwTT 2 (21) 

where wf p is the experimental window function, and n Q is a quantity denoting the noise 
level in ST/T per pixel. We compute Cf ET by combining the 45, 60 and 90 GHz channels 
of the MAP experiment. 

Next, we develop an optimal filter function which maximizes the signal-to- noise of the 
cross-correlation. 



3.1.1. Maximizing the signal 

We want to find the appropriate weight function, we(k,i]), that maximizes the signal 
from aimdim where, 

aim = afm+J H J k 2 dkA lm (k)j e {kr)')f{r,') (22) 
with the first term in the sum being the contribution from the surface of last scatter, and 



aim = I drj J k 2 dkA lm (k)D(r - T])w e (k,r])j e (kr]) (23) 
with Ai m (k) as the spherical harmonic expansion of the comoving number density, 

m{x) = I k 2 dkA lm (k)D(T - r])j e (kr])Y lm (n) (24) 



and 

m e c^ a(T — T]'Y 
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The noise term is 

([A lm (k)D(r - V )][A* lm (k)D(T - rf)]) = b'PWDiTo-vmTo-v^+nr^Wn^ir)') (26) 

where ^(77) depends on the galaxy selection function. In evaluating this function, space is 
considered split into a number of discrete time regions. We want to maximize 

"Kim 



( a Li)( a Z 2 m) 

[Ek k 2 J J drjdrj'wejk, v)ji(kr j )j e (kr l ')f(r l ')^P(k)D(r l )] 2 

(aL)Ek/ / dvdr]'w e (k,T])w £ (k,T]%(kr])j i (kT] , )[b^P(k)D(T])D(T)') + 

(27) 

Differentiating lnx z 2 m with respect to We(k,rj) and setting to zero yields: 

Ek k 2 J J d V d V 'w E (k, rj)jt{krj)j e (krf)f(rf)biP(k)D{rj) 

k 2 j £ (kT])[b 2 P{k)D{ri) j dri'jtjkri^DWwefari') + n^jjj) J d V 'j e {k V ')w £ {k, V ')] 

Ek II dr]dr]'w e {k, r])w e {k, r]')j e (kr])jfXkr]')[b 2 P(k)D(r])D(r]') + nj 1 ' 2 (v)n>i Xv')} 

(28) 

If we ignore the slow variations in the counts and growth factors, then the solution to this 
equation is 

* b 2 P(k)f( V >) 

We[ *> 71 } b 2 P(k)D(r - V >) + n-\ V ')/D(r - rf) ^ 

This weighting is just the projected signal over the projected noise. By weighting the 
integrand of Eq. [T^ with this optimal filter, one of the Cf i ~ sz factors in the numerator of 
Eq. cancels with the (C^ l ~ Ei + Cf olssON ) factor in the denominator, and we are left 
with: 

q v 2 A Q Q^i-SZ 

^ s r^V 2e+l) crrc^- (30) 



Prior to this point, the linear power spectrum P(k) was used in order to make the 
derivations of the cross-correlation and the optimal filter clearer. However, in the actual 
calculations we will use the fully-evolved non-linear power spectrum. In the non-linear case, 
the growth factor D gets absorbed into the power spectrum, which becomes a function of 
both wavenumber k and scale factor a: P(k, a). 

Let us define the quantity 

p(k,a) = y/P(k,a). (31) 
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Using this non-linear power spectrum, the quantities in Eq. RO are given by, 



% _ sz _ 2 (2a T k B n°Xb 



C, 



gas 



7T 



m P c z 



k 2 dk W s z (k)wf l - sz {k), 



with 



W f' sz {k) 
and Eq. O becomes 



o a(r - Tf) 
wf z (k) 



1 + b]P(k, a)ni{rj) 



o a(r - rj'y 



(32) 



(33) 



(34) 



Note that in the case of the SZ-ellipticity correlation for the FIRST objects, the 
signal-to-noise is defined analogous to Eq. O with 



a 



e-SZ 



2 (2a T k B n Q e T% 



6 6 u gas 



7i \ m e c 2 



k' 2 dkwf z {k)w e e (k) 



where 



Also, 



and 



w\{k) = / dr] 



TO 



1.68 x ^Qn ( 



#07) 



CI 



7T 



c / a(r - 77) 



a)jt{krj) 



a 



POISSON 



(35) 

(36) 

(37) 
(38) 



Nfirst 

where cr e = 0.4 as before, and N FIRS t is the number of FIRST objects per steradian. We 
estimate ~10 FIRST sources per square degree will be optically matched with the SDSS. 



3.2. ISW Effect 



The temperature fluctuations arising from the ISW effect can be expressed as, 



AT 



ISW 



drj 



<9$ 
drj 



We can relate the gravitational potential $ to the co moving density field 5: 

V 2 $(x) = ^5(x) 



(39) 



(40) 



- 12 - 



where the Laplacian is evaluated in comoving coordinates. Expanding 5(x) using Eq. 0, we 
obtain 



D(t -t)') r d 3 k 5(k) 



(r - rf) J (2vr) 3 / 2 k 2 



Combining Eqs. [39] and |4], the ISW term becomes 



AT 



m 



ISW 



3H$n r d 

' dr >dj 



a(r - rf) 



d 3 * S(k)^ 

(2tt) 3 / 2 fc 2 



kmr/ 



(41) 



(42) 



where we have used 8irGpo = 3HqQ . Following the same formalism as in § |3.1.1| to derive 
an optimal filter for the ISW-number-count cross-correlation, we obtain 



a 



7T 



f 2 wl sw {k) W f- ISW {k) 



where 



and 



w 



ISW 



w 



Xi-ISW 



(k) 



to Q 
or] 



KM 

a(r - rj) 



p(k,a') 
a(r - rf) 



je(kri') 



je{kr)) 



1 - 



l + blP(k } a)rii(ri) 



(43) 

(44) 
(45) 



Hence, following the notation of Eq. [TP], the signal-to-noise is: 

it,i-ISW 



ISW 



An 

Air 



C 



(46) 



We have neglected Cf ET in the denominator since the ISW effect is confined to low I where 
MAP's detector noise is very low. 



For FIRST objects, the signal-to-noise is 



Xi-ISW\2 



jsw 4vr V Cj T (Cf e + C[ OISSON ) 



where 



-isw 2 / 3i^Q VIq 



71 



dkwl sw {k)w\{k) 



(47) 



(48) 



and the other Cg and we terms are as defined previously in Eqs. |3B], |37|, and [0, 
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3.3. Weak Lensing 

We follow the methodology and notation of [Zaldarriaga fc Seljak (1999 )| . Most tracers 
of the underlying density field will not be perfect since they correlate only partially with the 
convergence k. Assuming we have a map Y that correlates with k and has a cross-correlation 
Cp Y , the cross-correlation between an e map and the Y map gives signal-to- noise 

I) 2 = *n y m,,(ef(cff 

N) M n, v (iw r .,(e) 2 cr + w)(c[ r + cf° ISSON ) k ; 

where e is a spin-zero quantity charactering the CMB polarization field. In this expression, 
the window function describing the effect of finite angular resolution is, 

ft d9 SjT ld£F(£)F(W + q- + 2£q cos 0)£ 2 Cj T 
WrM ~ 2nJ^£d£F(£) 2 Cr ( j 

where Cj T is the CMB power spectrum, and 

F(£) = 



i e<£cut 

£>£ C ut 



For the MAP experiment, £ cut = 600. Also, the large-scale amplitude of the noise spectrum 
is, 

2tt / fr* £d£tF{£)\Cj T + Cf ^ 
(/ £d££ 2 F{£) 2 Cj T y 
The convergence autocorrelation is given by 

2 



7T \ 2C 2 

where 

10? = 



c „ = |(3^«o\ (52) 



^ P(MI W)j<(fe/) (53) 
o a(r - 77') 



If the Y map consists of ellipticities of FIRST objects, then Cj Y = C| e , and C7 r becomes 



Cf = 2 |<3H|V ) / edk<(k)w e {k) (54) 



7T 



where w\ is from Eq. [36]. 



If the Y map consists of a number counts of tracer objects i, then we can use the 
optimal filter technique explained previously. Thus, the signal-to-noise becomes, 

\2 



S \ 2 AQ 4W re 



C, 



irS^ + fa 7^ . ^ (55) 



N J lensing 4vT ^ AW res (£) Cf K + N' ( 
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where 



and 



2 / 3Hq£Iq 



7T \ cr 



k 2 dkw K e (k)wf- K (k) 



.Si K 



, p k,a , , 
a(r - 7]') 



l + blP(k,a')ni(r)>) 



(56) 
(57) 



3.4. Cross-correlation Computation 

In evaluating the above integrals, the non-linear power spectrum P(k, a) was computed 
using the analytic approximation for a fully-evolved mass power spectrum for spatially flat 
cold dark matter cosmological models with quintessence given by [Ma et al. (1999)] . The 
linear power spectrum used in this approximation was computed using CMBFAST ( fSeljakj 
|fc Zaldarriaga 1996| ), using adiabatic initial conditions and a scale-invariant primordial 
spectrum. Going to a non-linear power spectrum makes very little difference to the 
cross-correlations, since most power arises from the linear regime, k <C 1 (h Mpc^ 1 ), at 
which the non-linear contribution to the spectrum starts to dominate. 

To increase the speed of the numerical integrations, the fc-integrals were done 
in logarithmic space and the spherical Bessel functions ji(krj) in the //-integrals were 
approximated by: 

. ( v / Hkrj) : kr < t 

3ei r])= \ yf%8(h,-£) : kr >i (58) 



The optimal filter technique was applied to galaxies (both the photo-z and complete 
samples) and to quasars with photometric redshifts. The number density functions ^(77) 
used in the number count optimal filters were computed as follows: 



1. The galaxy number density was calculated using the r-band Schecter-form luminosity 
function (LF) of galaxies in the Las Campanas Redshift Survey (LCRS) ([Lin et al.| 
1996Q . 

2. For photometric (i 1 < 21) quasars, the quasar luminosity function described in [Fau| 
(1999)| was used. This LF is defined for h = 0.5, where h = ifo/100 km s _1 Mpc -1 . 
Dependence on h was added to the LF using the fact that /i _3 $(M, z) and M — 5 log h 
are invariant under transformation of h. The LF was taken to be valid for absolute 
magnitudes -30.0 < M - 5 log(/t/0.5) < -23.0. Finally, this LF is defined in the B 
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band. To convert apparent magnitudes in the SDSS i' band to Mb, we have followed 
Schneider, Schmidt & Gunn (19911 



M B = i'-5 log [77(1 + z)} - 25.0 + 2.5 log(l + z) - 2.5a log 



7628 



;i + 2)4400 



0.15 (59) 



Here, rj(l + z) is the luminosity distance (in Mpc); a = —0.5 is the power law index 
of the quasar spectrum; the i' and B bands have effective wavelengths of 7628 and 
4400 A respectively, and the offset between the B and the AB magnitude systems is 
0.15 ( Fukugita et al. 1996; ). 



The comoving number density is given by, 



n(z) 



:i+z) 



where $(M, z) is the LF. 



$(M, z)dM 



(60) 



The cross-correlations of spectrometric quasars and FIRST objects were calculated 
using the non-filtered formalism (Eqs. |7 and IS). The radial distribution W(r]) of 
spectrometric {i 1 < 19) quasars was calculated using the redshifts of a set of quasars 
simulated in SDSS colors supplied by Donald Schneider. The radial distribution of FIRST 
objects was approximated by a Gaussian with mean z = 1.2 and width Az = 0.4. 



4. CROSS-CORRELATION RESULTS 

Table |I] shows the estimated (S/N) 2 that would be obtained by cross-correlating MAP 
with the SDSS survey for a pair of models with constant Q h 2 . 

The SZ and weak lensing cross-correlations give interesting values of signal-to-noise 
and can be used in various ways, described below. 



4.1. Tracing Gas Pressure Fluctuations 

While clusters produce the highest amplitude Sunyaev-Zeldovich fluctuations, large 
scale variations in density associated with the superclusters are likely the dominant 
source of large angle SZ fluctuations ( Persi et al. 1995| ; |Refregier et al. 1999| ). We can 
use the SDSS photo-z galaxy sample (r' < 21) to trace the large-scale structure out to 
z ~ 0.5. By measuring the cross-correlation between the MAP CMB measurements and 
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the galaxies in each photo-z slice, we can measure the amplitude of the gas pressure/galaxy 
cross-correlation in each photo-z redshift bin. 

For the SZ cross-correlation, we have from Eq. ^(J that (S/N) 2 oc b gas T®. The gas is 
highly biased {b gas ~ 4, |Refregier et al. 1999Q , and is at high temperature (T e ~ 0.5 keV). 
Since the surface number density of the photo-z galaxies is very high, the Poisson noise in 
the cross-correlation is minimized. Furthermore, the SZ signal is weighted by the growth 
factor but not by the scale factor for /3 — (Eq. |33]). Therefore most of the signal arises 
from low-redshift space where the galaxy density is highest. These factors contribute to a 
high signal-to-noise from this cross-correlation. 

Fig. [I] shows our predicted signal and error in our measurements of b gas p gas , where b gas 
is the ratio of the pressure fluctuations to the density fluctuations and p gas is the mean 
gas pressure. The points for b gas p are simply 1.936 x 10~ 16 a" 13 ergs cm -3 , which is just a 
constant when (3 = as we have here. This calculation illustrates how well the deviations 
from the average pressure can be constrained. The error bars are calculated from the inverse 
of the SNR, which is obtained from Eq. [T^, but with the tracer galaxies selected from a 
redshift " slice" . This is done by multiplying the radial distribution function by a gaussian 
of the required slice width, centered on the average redshift of that slice. Since most of the 
power in the cross-correlation arises in the linear regime, the covariance between redshift 
bins should not have a significant effect. 

These observations will yield the thermal history of the universe and a detection of the 
"missing baryons" ( pen fc Ostriker 1999| , L. A. Phillips et al. 2000, in preparation). 



4.2. Lensing and the CMB 

The weak-lensing signal is inherently small, so it is important to choose the low-z 
tracer carefully in order to maximize the SNR. Many different effects are important in 
determining how well a given class of low-z tracers will perform when cross-correlated with 
weak lensing of the CMB. 

Quasars are useful tracers of large scale density field. They are thought to be biased 
tracers of the underlying mass distribution. Since the quasar number counts are steep 
(din N/ dm ~ 1) for i' < 19, gravitational lensing by the foreground matter should produce 
significant fluctuations in the observed quasar angular correlations. These two different 
effects (bias and steep number-count slope) are important on different angular scales. 
Quasars can also probe the density field out to high redshifts (z ~ 4), thus sampling a 
significant fraction of the volume of redshift space. However, their surface number density 
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is low, so depending on the cut used to select the quasars, the Poisson error can be large. 

The surface number density of galaxies is very high, so the Poisson noise of the sample 
is small. However, the slope of the galaxy number counts is shallow (d In N/dm ~ 0.5 in r' 
band) so the intrinsic fluctuation term of Eq. ^ dominates over the lensing term and the 
cross-correlation with weak lensing does not trace the foreground matter fluctuations very 
well. The galaxy sample only extends out to z ~ 0.5, and thus probes only a small fraction 
of the volume probed by quasars. 

We find that cross-correlation of weak lensing with %' < 19 quasars in the SDSS survey 
is Poisson-noise-limited and does not provide very good signal-to- noise. The r' < 21 galaxy 
sample performs better, but is hampered by the fact that it does not extend out to large 
redshifts. Filtering i' < 21 photometric quasars, which have surface number density more 
than six times that of i' < 19 quasars (Fan 1999), with the optimal filter developed in 
§ [3.1.1| gives reasonably high SNR. We will discuss how to use this cross-correlation to 
constrain fl and the equation of state in § ||. 



4.3. Cross-correlations with Poor SNR 

The ISW effect does not produce a good SNR, because it dominates at low £ where the 
matter power spectrum is falling and also because the SDSS redshift survey does not probe 
to sufficient depth. Optically matched FIRST objects in the SDSS survey are Poisson-noise 
limited and are found to be useless for the purpose of cross-correlating with the CMB. 



5. EXPLORING THE EQUATION OF STATE 

The general equation of state of a density field p is 

p = poa- 3(1+u,) (61) 

where w is the ratio of pressure to density. For example, w = for dark matter, w = +1/3 
for radiation, and w — — 1 for A. An alternative to the cosmo logical constant A as the 
missing energy in a flat universe with Q < 1 is given by quintessence Q ( |Caldwell, Dave 



fe Steinhardt 1998|) . It is a time-evolving, spatially inhomogeneous energy component with 



negative pressure and an equation of state w < 0. It drives cosmo logical expansion at late 
times, influencing the growth of structure. On very large scales, the quintessence clusters 
gravitationally, thereby modifying the level of CMB temperature anisotropy relative to 
the matter power spectrum amplitude and also inducing a late-time ISW effect. On small 
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length scales, fluctuations in Q disperse relativistically and the Q-field behaves as a smooth 
component. Quintessence modifies the matter power spectrum, the time evolution of the 
scale factor a and the linear growth factor D. In this section we investigate the possibility 
of detecting observational imprints of Q by cross-correlating MAP and SDSS. 

Following the notation of § |3.3| , the x 2 f° r cross-correlating weakly- lensed i' < 21 
quasars with the CMB convergence is given by, 



K (n ,w,b q )-cf<-"(w ,w',v q ) 



C * K + N ee/ m ^f\ Cf q - K {W , W>, b' 



c 



ffjPOissoNy 



(62) 



We take the fiducial models Cf q K (Q' ,w' \b' 



and Cf"~* q (n' ,w',b' q ] 



to have 



[Qo = 0.35, h 2 = 0.7,w = —1.0,b q = 3.0]. It is seen from Figure ^ that there is a degeneracy 
between the quasar bias b q and Qq. For each model, we minimize \ 2 with respect to b q . 

Exploring the parameter space for a family of spatially flat models with constant Q h 2 
for weak lensing of i! < 21 quasar number counts, we obtain the likelihood surface shown in 
Figure |3|. Here, we have assumed that there are 60 i' < 21 quasars/sq. deg., in accordance 
with [Fan ( 199971 - T ne labeled curves project to one-dimensional intervals containing 68.3%, 
Vo and 99% of normally-distributed data. 



At 90% probability, we obtain the following bounds: 0.30 < Qo < 0.42 and 
-1.0 < w < -0.77. 



6. SUMMARY 

We apply a formalism for computing the cross-correlation between data-sets at high- 
and low-redshift to MAP data and various classes of tracer objects from the SDSS, for 
a range of spatially-flat low-f2o models with a general equation of state. Three low-z 
effects that imprint secondary anisotropics on the CMB are considered: the thermal 
Sunyaev-Zel'dovich effect, the Integrated Sachs- Wolfe effect and weak lensing. 

Using photometric redshifts, one can select a tracer sample to match the redshift 
windows where a particular second-order CMB anisotropy is dominantly formed, using the 
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optimal filtering technique detailed here. As a practical algorithm, one would derive the 
luminosity function of the tracer sample, and hence the comoving number density of that 
sample for a particular cosmology. Then one would guess a power spectrum and tracer 
bias, making use of the CMB and tracer power spectra; construct the optimal weighting for 
each 3-D Legendre mode as a function of redshift and multipole; and finally project to form 
a two-dimensional map which is then cross-correlated with the CMB. One then estimates 
the power spectrum and bias from this map, and repeats the procedure until the process 
converges to the real optimal filter. 

We find that these techniques would give good constraints on pressure fluctuations 
of supercluster-scale gas and provide constraints on Q and the equation of state w in 
quintessence models. It is found that the photometric-z quasar and galaxy samples expected 
from the SDSS would perform well as low-z tracers. However, the complete SDSS survey 
would be required to generate sufficient signal-to-noise. We also find that the ISW effect 
does not provide enough signal-to-noise to make use of these methods. 

We are very grateful to Donald Schneider for providing the radial distribution of 
quasars with %' < 19 from his SDSS simulation; Alexandre Refregier for useful discussions 
on the SZ effect, cross-correlating FIRST ellipticities, and the bias of the intracluster gas; 
Xiaohui Fan and Gordon Richards for information on the SDSS quasar sample; Michael 
Strauss, Daniel Eisenstein, Andrew Connolly and Rita Kim for information on the SDSS 
galaxy sample; and Mark Jackson for proofreading this manuscript. HP and DNS were 
partially supported by the MAP/MIDEX program and NASA NAG5-7154. 
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Table 1. (S/N) 2 from the Cross-Correlations. 



(S/N) 2 (h = 0.7, Vl = 0.35, A = 0.65, T° = 0.5 keV) 


SZ Effect 


Weak Lensing 


ISW Effect 


Photometric galaxies (r 1 < 21, b g = 1) 


140 


29 


6 


All galaxies (r' < 23, b g = 1) 


175 


77 


10 


Spectrometric quasars {i' < 19, b q = 3) 


4 


17 


3 


Photometric quasars (i f < 21, 6 g = 3) 


31 


125 


8 


FIRST objects (6/ = 2) 


< 1 


< 1 


< 1 


2MASS galaxies a 


108 


3 


1 


(S/A0 2 (h = 0.76,^0 = 0.3, A 


= 0.7,T e ° = 


0.5 keV) 




Photometric galaxies (r' < 21, b g = 1) 


124 


25 


7 


All galaxies (V < 23, 6 9 = 1) 


158 


67 


13 


Spectrometric quasars {i' < 19, b q = 3) 


4 


15 


3 


Photometric quasars (i r < 21, 6 g = 3) 


31 


119 


11 


FIRST objects (6/ = 2) 


< 1 


< 1 


< 1 


2MASS galaxies a 


93 


2 


1 



a Estimated signal-to-noise values for the 2MASS survey flChester et al. 1998| ) are given as a 
comparison. This all-sky survey will contain < 10 6 galaxies brighter than K s ~ 13.5 (corresponding 
roughly to r' ~ 17), and about three quarters of the sky will be useful for the purpose of cross- 
correlating with the CMB. The complete 2MASS dataset will be available much sooner than that 
of the SDSS. 
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Fig. 1. — Pressure Fluctuations This figure shows the predicted measurement of 
the comoving pressure fluctuations based on measurements of the correlation between 
photometrically-selected galaxies (r' < 21), separated into redshift "slices", and CMB 
temperature fluctuations (see § |4.1|) . Pressure is in cgs units. The calculation assumes 
that T e = 0.5 keV today and b gas = 4. 
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Fig. 2. — Quasar-Lensing Cross-Correlation This figure shows the predicted correlation 
between %' < 21 quasars in the SDSS survey and the convergence computed from the CMB 
sky (see § O). The solid lines are for a ACDM model with fi = 0.35, A = 0.65, H = 70, 



and b q = 2. The error bars are based on two year MAP data and the full SDSS survey. The 
dotted lines are a model with Q = 0.3, A = 0.7, H = 76 and b q = 2. The dashed lines are 
for the same model with b q = 3. The degeneracy of the parameters fl and b q is apparent. 



-25 - 



a 0.40 



0.60 



0.50 




0.30 



0.20 



--1.00 -0.90 -0.80 -0.70 -0.60 -0.50 

w 

Fig. 3. — Parameter Estimation This figure shows the x 2 surface for the f2 — w parameter 
space of a family of spatially flat models with constant £l h 2 . These constraints are due to 
cross-correlation of weak lensing of quasar number counts for i' < 21 with MAP CMB data. 
For each model, \ 2 has been minimized with respect to the quasar bias b q . The fiducial 
model [O = 0.35, h 2 = 0.70, w = —1.0,b q = 3.0] is indicated by the diamond. The curves 
project to one-dimensional intervals containing the probabilities with which they are labeled. 



